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ABSTRACT 
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_C Aims. We present INTEGRAL light curves and spectra of the black-hole binary Cyg X-1 during a bright event that occurred in 2006 

0- September, and which was simultaneous with a detection at 0.15-1 TeV energies by the MAGIC telescope. 

Methods. We analyse the hard X-ray emission from 18 to 700 keV with the INTEGRAL data taken on 2006 September 24-26 by the 
IBIS and SPI instruments. These data are supplemented with RXTE All Sky Monitor data at lower energy. We present the light curves 
and fit the high energy spectrum with various spectral models. 
Results. Despite variations in the flux by a factor of ~2 in the the 20-700 keV energy band, the shape of the energy spectrum remained 
remarkably stable. It is very well represented by an e-folded power law with the photon index of T = 1 .4 and a high energy cut-off 
at £ c a 130-140 keV. The spectrum is also well described by thermal Comptonisation including a moderate reflection component, 
, with the solid angle of the reflector of ~ 0.4 x 27r. The temperature of the hot Comptonising electrons is kT e ~ 70 keV and their 

^ • Thomson optical depth is t ~ 2.5. These spectral properties are typical of those observed in the low/hard state. This shows that Cyg 

X- 1 may stay in the low hard state at least up to the flux level of 2 Crab, which corresponds to ~2-3 % of the Eddington luminosity. 
It is the first time a persistent high-mass black-hole binary is observed at a few percent of the Eddington luminosity with a stable 
low/hard state spectrum over a period of a few days. Such a bright hard state has so far been observed only during the rising phase 
■ of transient low-mass black-hole binaries. The TeV detection coincides with the peak of a small X-ray flare just after a very fast rise 

in hard X-ray flux. In contrast, the source remained undetected by MAGIC at the peak of a larger X-ray flare occurring one day later 
and corresponding to the maximum of the X-ray luminosity of the whole outburst. We do not find any obvious correlation between 
the X-ray and TeV emission. 
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?H \ 1. Introduction dred keV or so (e.g., Gierliriski et al. 1997). Occasionally, the 

. . , source switches to the high/soft state. The high-energy power 

Cyg X-1 is the prototype of black hole binaries. Since its dis- law is then much softer (r > 2 4) and the bolomet ric luminos- 

covery in 1964 (Bowyer et al. 1965), it has been intensively ob- ity is dominated by a thermal component peaking at a few keV 

served by all the high energy instruments, from soft X-rays to y- (e g #> Gierliriski et al. 1999). In recent years, Cyg X-1 also of- 

rays. It is a persistent source, powered by accretion onto a black ten ap peared in intermediate states in which the source exhibits 

hole from a massive companion, HDE 226868, most likely via a spectra that are interme diate between the two canonical spectral 

focused wind. The value of the mass of the black hole is subject states (re l a tively soft hard X-ray spectrum with T ~ 2.0-2.3) and 

to controversy, it is in the range of M x - (5-15)M© accord- with st flari activit Malzac et al (2006) 
ing to Herrero et al. (1995) or M x * (14-27)M© (Ziolkowski 

2005; see also Gies & Bolton 1986). The distance, D, is most INTEGRAL is regularly observing Cyg X-1 during monitor- 
likely within D 2. 1 + 0.2 kpc (Ziolkowski 2005 and references in 8 observations of the Galactic Plane. On 2006 September 24- 
therein) - hereafter we adopt D - 2 kpc 26, the source was at its highest hard X-ray level since the launch 
This source is most often observed in the so-called low/hard of INTEGRAL in 2002 October (Tiirler et al. 2006). During the 
state, characterised by a relatively low flux in the soft X-rays (~1 two " da y observation, the flux of the source varied significantly 
keV) and a high flux in the hard X-rays (-100 keV). In the hard and reached values > 2 Crab ' Further INTEGRAL observations 
state, the high energy spectrum can be roughly described by a of C yg X " 1 taken on 2006 September 30 indicated that the source 
power-law with photon spectral index, T, varying in the range returned to a more usual state at a level of ~1 Crab in the 20-40 
~ 1 .4-2, and a cut-off at a characteristic energy, E B , above a hun- keV band and 1 3 Crab ln ^ 4 °- 80 keV band - 

The peak fluxes measured during 2006 September 24-26 

Send offprint requests to: J. Malzac were among the highest daily fluxes observed by the CGRO 
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Fig. 1. Long term hard X-ray light curves of Cyg X-l. The 
SWT/BAT (15-50 keV) and the RXTE/ASM (1.5-12 keV) 
light curves are displayed in panel (a) and (b) respectively. The 
BAT Crab unit is 0.221 crrT 2 s _1 , and the conversion of the 
ASM counts into physical flux was done following Z02. Panel 
(c) shows the evolution of the photon spectral index based on the 
3-5 and 5-12 keV bands, as described in Z02. The red arrows at 
MJD 54003 point at the time of the INTEGRAL observation de- 
scribed in this paper. The blue arrows at MJD 54604 indicate the 
event reported by Neronov et al. (2008), and the green arrows 
show the recent short flare peaking at MJD 56900. 



Burst and Transient Source Experiment (BATSE) from 1991 to 
2000 in the 100-300 keV band and exceeded all of them in the 
20-100 keV band (Zdziarski et al. 2002, hereafter Z02). Similar 
high fluxes were also visible in the SWIFT Burst Alert Telescope 
(BAT) and RXTE All Sky Monitor (ASM) light curves. Fig. [TJ 
shows the long term evolution of the daily fluxes observed by the 
BAlQ and ASM. The event detected by INTEGRAL coincides 
with the highest flux ever measured by BAT (15-50 keV) since 
the launch of SWIFT in 2004. The evolution of the X-ray spectral 
index, shown in Fig. [Tic), demonstrates that Cyg X-l was then 
in the hard spectral state, which started around MJD 53900, and 
has lasted since then till now (MJD 54720). The ASM (1.5-12 
keV) fluxes were also among the brightest ones in the hard state. 
We refer to Albert et al. (2007, hereafter A07) for additional 
close-ups of the BAT and ASM light curves around the time of 
the INTEGRAL event. Based on data from the ASM, the X-ray 
flareQ lasted about three days, about MJD 54002.0-54005.0 (see 
figure 4 of A07). 

Even brighter flares have already been reported (Stern et 
al. 2000; Golenetskii et al. 2003; Gierlinski & Zdziarski 2003), 



1 http://swift.gsfc.nasa.gov/docs/swift/results/transients/CygX-l 

2 A flare is a change in intensity of the source with, typically, a ris- 
ing phase followed by a decay phase. It is convenient and customary to 
describe the X-ray variability of X-ray binaries in terms of a superposi- 
tion of flares of various amplitudes and time scales. We stress that this 
description does not presuppose any physical mechanism for the origin 
of this variability, nor that each of these flares is genuinely related to an 
individual localised physical event. 



however these events were all of much shorter duration (ms to 
ks rather than days). The 2006 event represents the first obser- 
vation of a persistent black hole binary in a stable hard state at 
such levels of intensity over such a long time. This event is also 
particularly interesting with regards to the simultaneous, and the 
first ever, detection of Cyg X-l at TeV energies by the MAGIC 
telescope on 2006 September 24 (MJD 54002.928-982, A07). 
The 0.1-1 TeV luminosity during the MAGIC observation was 
^ 2.2 x 10 34 erg s _1 , whereas the X-ray flare luminosity was 
about 5 x 10 37 erg s _I . The TeV detection was during the or- 
bital phases of 0.90-0.91 (using the ephemeris of LaSala et al. 
1998; Brocksopp et al. 1999), with the black hole behind the 
companion. Also, as pointed out by Poutanen et al. (2008) and 
Rico (2008), Cyg X-l was then close to the the peak of the long- 
term, superorbital, flux modulation (seen both in radio and X- 
rays, e.g., Lachowicz et al. 2006). 

This paper presents a detailed analysis of the INTEGRAL 
light curves and spectra obtained during this bright event. 



2. Data analysis 

We have analysed the data corresponding to the rev. 482 (MJD 
54002.4-54004.3). During this revolution, INTEGRAL was per- 
forming a Galactic Plane scan observation, and for this reason 
Cyg X-l was always seen at least 6° from the centre of field 
of view. Therefore, the data are not optimal. The effective ex- 
posure was decreased by a factor of ~3 when compared to the 
standard dithering observation of the source. As Cyg X-l was 
always outside the JEM-X field of view and no other X-ray tele- 
scope pointed at the source during this event, we lack spectral in- 
formation at S: 18 keV, except for the ASM fluxes in three bands 
within 1.5-12 keV, see Section[3] In total, there were 51 point- 
ings in rev. 482 when Cyg X-l was observed at < 18° off-axis. 
However, to avoid systematic errors appearing for large off-axis 
observations we finally chose only 3 1 pointings where the source 
illuminated at least 10% of the ISGRI detector surface (with the 
off-axis angle between 6° and 13°). 

The IBIS/ISGRI and SPI data were reduced with the stan- 
dard software of the Off-Line Scientific Analysis package (OSA, 
v. 7.0, 2007 September) released by the INTEGRAL Science 
Data Centre (ISDC, Courvoisier et al. 2003). For ISGRI data 
analysis, we used the default OSA 7.0 settings and default re- 
sponse files. In case of SPI the default flat field option was cho- 
sen for the background modelling, with the constant background 
interval set to 5 pointings. For both instruments, the input source 
catalog contained all strong objects in the Cyg X-l field, namely 
Cyg X-l, KS 1947+300, 4U 1957+1 15, QSO B1957+405, EXO 
2030+375, Cyg X-3, SAX J2103.5+4545. The selection was 
based on the ISGRI mosaic image. 

IBIS/PICsIT data analysis was performed with the non- 
standard tools developed for count rate extraction from spectral- 
imaging, SINGLE events data (Lubinski 2008). Basic elements 
of that method are: direct handling of Poisson probability density 
functions for both background map normalization and source 
count rate, careful modelling of the background variability and 
energy-dependent pixel illumination functions. This new ap- 
proach allowed for a detection of several weaker sources and 
improved the results for stronger sources when compared to the 
standard OSA PICsIT software results (Lubinski 2008). Usually 
PICsIT background maps are made for entire revolution but 
when some background evolution is observed, as in the case of 
rev. 482, a better result is obtained when there are more back- 
ground maps. Two maps, prepared for periods before and after 
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Fig. 2. IBIS/ISGRI light curves during rev. 482 (MJD 54002.4-54004.3) The upper panel shows the light curves in the 18-40 keV 
(triangles), 40-100 keV (squares) and 100-200 keV (stars) bands. Each point represents the time averaged flux over a science 
window (a pointing of duration of ~2 ks). For clarity, the data points of the 18-40 keV, 40-100 keV and 100-200 keV light curves 
are connected by dotted, dashed and dash-dotted lines, respectively. The three unconnected symbols (triangle, square and star), 
around 1 Crab, show the historical average IBIS/ISGRI fluxes of Cyg X-l in the respective energy bands. The average fluxes were 
calculated using all the pointings with Cyg X-l in the field of view since the launch of INTEGRAL. The long vertical lines associated 
with these symbols show the amplitude of variability (rms) of the historical INTEGRAL light curves. The bottom panel shows the 
hardness ratios of the count rate in the 40-100 keV band to that in the 18-40 keV band (squares), and the corresponding 100-200 
keV to 20^1-0 keV ratio (stars). The dark vertical stripes indicate the times while the source was observed with the MAGIC telescope. 
The thick vertical dashed lines mark the start and the end of the MAGIC pointing in which Cyg X-l was detected at TeV energies, 
the other pointings showed no detection. The two vertical solid lines enclose the two science windows corresponding to the peak 
spectrum shown in Fig. [5] 



the end of science window number 3 1 , allowed us to separate the 
background variability from the source variability. 

Hereafter, the light curves are presented in Crab units (see 
Figs. [2]-|4]i. For IBIS/ISGRI, we used the standard count rates 
provided by OSA and then define 1 Crab as 179.4, 98.3 and 17.7 
s 1 in the 18^10, 40-100 and 100-200 keV energy bands re- 
spectively. For the ASM, we used a Crab count rate of 75 s _1 . 
For the PICsIT, 1 Crab was defined as 3.44 s _1 in the 277-332 
keV band. The PICsIT light curve, shown in Fig. [3] is probably 
the first hour-scale light curve for any source prepared for a sim- 
ilar energy band. Cyg X-l is detected in practically all science 



windows because the corresponding 3cr upper noise limit for 3.5 
ks observation, estimated on a basis of Monte Carlo simulations, 
is about 0.5 Crab. 

We have produced three spectra for each of the three instru- 
ments (IBIS/ISGRI, PICsIT and SPI). The first one is the total 
spectrum, i.e., averaged over the entire revolution, but with the 
pointings only for the offset < 13° (see above). We used the same 
science windows (pointings) for all instruments, which resulted 
in 31 science windows distributed over the rev. 482 with total 
exposure =s 70 ks. The second spectrum is the average of the two 
pointings around the peak of the flare, (pointings 46 and 47) with 
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the exposure of 5.2 ks, see Fig. [2] The third one is an average 
of the two pointings spanning the time period MJD 54002.924- 
54003.007 coincident with the MAGIC detection at TeV ener- 
gies. 
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Fig. 3. Light curves of Cyg X-l during rev. 482. The upper pan- 
els show the IBIS/ISGRI light curves in the 18-40, 40-100 and 
100-200 keV bands with a resolution of 100 s. The two bottom 
panels show the IBIS/PICsIT light curve in the 277-362 keV 
band, and the RXTE/ASM 1.5-12 keV light curve. 

The spectral analysis was performed under XSPEC v. 11.3 
(Arnaud et al. 1996). The energy range used was 18-300 keV 
for the IBIS , 24-700 keV for the SPI, and 277-632 keV for the 
PICsIT. A 1.5% systematic error was added in quadrature to the 
statistical errors of the IBIS/ISGRI spectra. Moreover, in order to 
account for uncertainties in the cross-calibration when simulta- 
neously fitting spectra from different instruments, a multiplica- 
tive constant was added in the spectral fits to each instrument 
data set: it was set free for IBIS and frozen to 1 for SPI. The 
resulting best-fit value for these constant was around 0.9 for 
ISGRI. For PICsIT the normalisation constant is not very well 
constrained by the data. We found values in the range 0.8-1.2 
depending on the data set and model used. 
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Fig. 4. 100-s resolution IBIS/ISGRI light curves of Cyg X-l 
around the time of the MAGIC detection. The triangles and 
squares stand for the light curves in the 18—40 keV and 40-100 
keV band respectively. The time span of the MAGIC detection 
is delimited by the two vertical dashed lines. 
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Fig. 5. SPI spectra at the outburst peak (blue) and averaged over 
the rev. 482 (red), fitted by the e-folding power-law model (see 
Table [TJ. For comparison, we also show the best-fit model of a 
typical hard-state spectrum (period 1 of Cadolle Bel et al. 2006). 



3. Results 

The IBIS/ISGRI light curves of rev. 482 are shown in Fig. [2] As 
shown on this figure, during the whole observation the measured 
fluxes were a factor of 2 above the historical average IBIS/ISGRI 



fluxes and also well outside the range of the average variability 
amplitude. The light curves can be decomposed into an under- 
lying approximately constant component at ^1.5 Crab (18-40 
keV), =d.8 Crab (40-100 keV) and =T.9 Crab (100-200 keV), 
and, on top of it, three main flares lasting for several hours and 
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Table 1. The best fit parameters of the Cyg X-l spectra when fitted with an e-folded power law; the e-folding energy, E c , the photon 
index, F, and the normalization at 1 keV, K, in units of photons keV~' cirT 2 s _1 . The results are shown for the spectra averaged over 
rev. 482 and for the peak of the outburst (as indicated). In addition, the last line shows the results for the spectrum averaged over the 
2 pointings that were simultaneous with the MAGIC TeV detection. 
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Fig. 6. Simultaneous INTEGRAL/MAGIC spectrum of Cyg X- 
1. The SPI (blue), ISGRI (cyan) and PICsIT (red) spectra are 
averaged over the time of the TeV detection. For clarity, all 
the INTEGRAL data were rebinned. The ISGRI data were nor- 
malised to the SPI spectrum (the ISGRI flux was divided by 
0.96) as determined from the fit with a cut-off powerlaw model 
(see TableHJ. The black crosses show the simultaneous MAGIC 
spectrum (A07). The dashed line is a power-law with photon in- 
dex F = 2.5 connecting the two data sets. 
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reaching fluxes > 2 Crab. The strongest outburst and reached the 
maximum of about 2 Crab (18-40 keV), 2.4 Crab (40-100 keV) 
and 2.5 Crab (1 00-200 keV) in a pointing on 2006 September 25 
20:58-21:57 UT. This main flare lasted for 10 hours (estimated 
from the time interval between the 2 minima in the light curve 
before and after the peak). 

The times when MAGIC was pointing at the source and 
the time of the TeV detection are also marked on this figure. 
Strikingly, MAGIC did not detect Cyg X-l at the peak of the 
10-hour outburst (A07), but one day earlier, when the the hard 
X-ray flux was close to the average level of our observation. 

The 40-100 keV to 18^10 keV (hereafter HR1) and the 
100-200 keV to 18-40 keV (hereafter HR2) hardness ratios are 
shown in the bottom panel of Fig. [2] Both hardness ratios are 
remarkably stable: the spectrum remains impressively constant 
in spite of the relatively strong flux fluctuations. Nevertheless, 
although we do not detect any fluctuation of HR1, we find small 
(<15%) but significant variations of HR2. As seen on Fig.|2j the 
HR2 tends to be anti-correlated with the flux. This is confirmed 
by a Spearman rank test showing that HR2 is anti-correlated 
with the flux in the 40-100 keV energy band with probability of 
2 x 10~ 6 that this correlation occurred by chance. This behaviour 
is similar to that shown during an INTEGRAL observation during 
the normal hard state, see fig. 3 in Bazzano et al. (2003), and, on 
longer time scales, by the BATSE, see fig. 6a in Z02. The origin 



Fig. 7. SPI (blue), ISGRI (red) and PICsIT (green) spectra aver- 
aged over rev. 482 fitted with the model consisting of thermal 
Comptonisation (dashed curve) and Compton reflection (dot- 
dashed curve) model, see Table [2] The solid curve gives the 
sum. The triple-dot-dashed curve shows the same model but 
taking into account the standard absorption column density to- 
ward the source (Nh = 6 x 10 21 cm -2 ). For clarity, all the data 
were rebinned and normalised to the SPI spectrum (the ISGRI 
and PICsIT fluxes were divided by 0.889 and 1.18 respectively). 
The three dark horizontal lines below 12 keV show the average 
ASM energy fluxes overMJD 54002.4-54004.3 (those data were 
not taken into account in the fit procedure). 



of this effect was modelled by Z02 as due variability of the total 
luminosity in a self-consistent thermal Comptonization model. It 
is also worth noting that HR2 was at a very similar level during 
the MAGIC detection and at the peak of the outburst. Thus the 
TeV emission seems unrelated to the hardness ratio. 

Fig.[3]shows the IBIS/ISGRI light curves with a higher time 
resolution, 100 s, which reveal a strong flaring hard X-ray ac- 
tivity also occuring on shorter time scales. Using the full light 
curve to compute the rms amplitude of the variability, we find 
that the Poisson noise corrected rms variability amplitude is of 
23, 27 and 23% in the 18^0, 40-100 and 100-200 keV bands 
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respectively. Statistical noise contributes at a level of 3.2, 5.5 and 
51%, respectively, to the variance of the plotted light curves. 

The intrinsic variability appears, however, dominated by the 
longer time scales. In order to evaluate the amplitude of vari- 
ability occurring on time scales shorter than one hour, we com- 
puted the rms on segments of duration 3400 s, which removes 
the longer time scales. The average, noise corrected, rms ampli- 
tude is then reduced to 7% in all three energy bands. The noise 
contributes at a level of 10, 16 and 73% to the observed variance 
(i.e. most of the short time scale variability of the 100-200 keV 
band seen on Fig. [3] is actually due to statistical noise). 

A close-up of the light curves around the time of the MAGIC 
detection is presented on Fig. [4] The start of the TeV emission 
coincided with the peak of an X-ray flare, while there was no 
TeV detection for 2.5 hours before. However the MAGIC detec- 
tion could also be seen (from Fig. [3]and|4]i as preceeding the 
sharp rise of a flare not fully covered by INTEGRAL observa- 
tions. If the TeV emission is indeed a precursor and not simul- 
taneous to the hard X-ray emission, this may help to understand 
the absence of TeV detection during the peak of the major 10- 
hour flare. Alternatively, we note that the detection occurs just 
after a rapid (< 5 min) increase in the 40-100 keV X-ray inten- 
sity. This brightening appears to be the most abrupt change in 
flux during the whole INTEGRAL observation (see Fig. [3} and 
the TeV emission might also be associated to this fast rise in the 
X-ray luminosity. In fact, the lack of coverage and the unique- 
ness of the TeV detection prevent us from linking firmly the TeV 
emission to any feature of the hard X-ray light curve. 

Fig- E] shows both the spectrum averaged over rev. 482 and 
the peak SPI spectrum. In addition it shows, as reference, the 
best-fit model of a typical hard-state spectrum previously ob- 
served by INTEGRAL (Cadolle et al. 2006). Notably, it has a 
flux lower by a factor of ~2 but a very similar spectral shape. 
Both SPI spectra are well represented by an e-folded power-law 
model characterised by a power-law photon index F ~ 1 .4 and a 
cut-off energy E c ~ 140 keV. The best fits parameters are shown 
in Table Q] Both spectra have the best fits parameters consistent 
with being the same within the uncertainties, confirming the ab- 
sence of any strong spectral evolution during the flare. We have 
also performed similar fits using the ISGRI spectra and com- 
bining the SPI, ISGRI and PICsIT. The best fit parameters are 
shown in Table Q] and are similar to those obtained using SPI 
alone. The small differences in the best fit parameters are due to 
known calibration problems. The spectrum during the outburst 
is consistent with a stable hard state. This shows that Cyg X-l 
may stay in the hard state at least up to fluxes of 2 Crab. 

Fig. [6] shows the spectrum of the source at the time of the 
MAGIC detection, together with the simultaneous TeV data 
(A07). Although the statistics is relatively poor, the INTEGRAL 
spectrum appears to be very similar to the average and peak 
spectra (see best fit parameters in Table[TJ. 

We also fit the average spectrum with the Comptonisation 
model compps (Poutanen & Svensson 1996). In this thermal 
Compton model, the spectrum is computed for a homogenous 
spherical hot corona, inside which blackbody photons are in- 
jected and then Compton upscatered by hot Maxwellian elec- 
trons. The corona is parametrized by its electron temperature, 
T e , and the radial Thomson optical depth, t = n e <x x //, where n e 
is the electron density, crj is the Thomson cross section, and H 
the radius of the sphere. The temperature of the seed photons is 
kept fixed at 0.3 keV. This model also accounts for the reflection 
of the hard X-ray radiation illuminating the disc and forming an 
additional reflection bump peaking at ~30 keV (Magdziarz and 
Zdziarski 1995), for which we assume the inclination of 45°. The 



Table 2. The best fit parameters of the Cyg X-l spectra for 
thermal Comptonization; the electron temperature, kT e , the 
Thomson optical depth, t, and the strength of Compton reflec- 
tion, R, relative to that from an infinite slab. The fit results with 
the fixed (denoted by 'f') R — are also shown. 



instrument 


kT e (keV) 


T 


R 




SPI 


64+ 3 


2 36 +u uv 

2 25 +im 
-0.09 
2 94+0.15 

2 61+^ 

? 7^+0.14 
A - -0.05 


Of 


2.51(22) 


SPI 
ISGRI 


75+' 
-S 
55i 


0.40 + °* 
Of ' 


1.24 (21) 
2.36(48) 


ISGRI 
SPI+ISGRI+PICsIT 


64i 
58+' 


26 +01 
Of 


1.14 (47) 
2.54(78) 


SPI+ISGRI+PICsIT 


68i 


2 46 +0 07 

■^"-0.07 


30 +009 


1.13 (77) 



strength of the reflection component is quantified using its rela- 
tive amplitude R with respect to the case of the X-ray emission 
of an optically-thin isotropic source reflected by an infinite slab 
of cold material. The abundances of the reflector are solar and 
the reflecting material is neutral. 

The best fit parameters are shown in Table [2] We fit the 
ISGRI and SPI spectra independently as well as the combined 
ISGRJ/PICsIT/SPI spectrum. In the fits, we find kT t ~ 70 keV 
and t ~ 2.5. The presence of a reflection component is required. 
The fits performed without it are poor, with^ 2 /v > 2.3, whereas 
its inclusion improves the fits dramatically, leading to^- 2 /v 1. 
The combined ISGRI/SPI/PICsIT spectrum, shown in Fig. [7] re- 
quires a moderate reflection strength of R ~ 0.4. Overall, the pa- 
rameters of the thermal Compton and reflection model are sim- 
ilar to those usually obtained in the low/hard state (see, e.g., 
Done et al. 2007). In Fig. [7] we also show the average ASM 
energy fluxes over MJD 54002.4-54004.3, converted from the 
count rates using the response matrix of Z02. These fluxes turn 
out to be consistent with the low energy extrapolation of the best 
fit model of the INTEGRAL spectrum, confirming the validity 
of the Comptonization model (the extrapolation of the best fit e- 
folded power law model does not provide such a good agreement 
with the ASM data). 

Non-thermal high energy components in excess of thermal 
Comptonisation models have been reported in the hard state of 
Cyg X-l (McConnell et al. 2000, 2002, Cadolle Bel et al. 2006), 
and, e.g., the low-mass X-ray binary GX 339^1 (Wardziriski et 
al. 2002; Joinet et al. 2007). Such a hard tail is usually attributed 
to the presence of a small fraction of non-thermal electrons in 
the hot comptonising medium. In the present data the presence 
of such an excess is suggested by the shape of the SPI resid- 
uals above 300 keV (see in Fig. |7). We therefore fit again the 
combined spectrum shown in Fig. [7] allowing for the presence 
of a non-thermal tail in the Comptonizing electron distribution. 
This did not lead to any sensible improvement of the x 1 - This 
excess does not appear in the PICsIT data, but this depends on 
the poorly constrained normalisation factor between SPI and 
PICsIT. We therefore conclude that, even though it is not ex- 
cluded, we do not find robust evidence for the presence of a non- 
thermal excess in this observation. 



4. Discussion 

Events similar to the luminous state presented in this paper were 
previously reported in Cyg X-l. Stern et al. (2001) discovered 
bright hard X-ray flares lasting ~1 ks in the BATSE data. These 
events occured in the hard state with the peak luminosity at > 30 
keV reaching 1.6 x 10 38 erg s _1 , which is an order of magni- 
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tude larger than the normal luminosity of Cyg X-l in this state. 
They found two types of spectral behaviour during these out- 
bursts, one with constant spectrum, similar to what we observed 
in September 2006, the other showing spectral varibility sug- 
gesting the presence of two spectral components. Golenetskii et 
al. (2003) reanalysed the archives of the BATSE, Ulysses, and 
the Wind/Komis experiments, and found 7 strong outbursts of 
duration < 28 ks, both in hard and soft state, with the peak 
15-300 keV luminosity similar to the one found by Stern et 
al. (2001). Fig. [8] shows the spectrum fitted to the brighest of 
those flares (in the hard state) and compares it with the peak 
INTEGRAL/SPl spectrum, as well as with the average Cyg X-l 
hard-state spectrum from McConnell et al. (2002). Note that the 
peak flux of brightest previous flare was still a factor 1.8 higher 
than that shown in Fig. [8] (see fig. 5 and compare tables 2 and 3 
in Golenetskii et al. 2003). The isotropic bolometric luminosity 
inferred from that fit and adjusted to the peak is L ~ 3 x 10 38 
erg s , which is ~ 2O%(1OM0/M) of the Eddington luminos- 
ity, Le — 1.5(M/10Mq)10 39 erg s _1 (assuming the solar abun- 
dances). 

In addition, powerful sub-second flares were found by 
Gierlihski & Zdziarski (2003; see also Zdziarski & Gierlihski 
2004), mostly in the hard state. The hard-state flares reached 
the 3-30 keV luminosity of ~ 10 38 erg s _1 . Their estimated 
bolometric luminosity was up to ~ 4 x 10 38 erg s 1 , which is 
~ 25%(10Mq/M)Le- We have searched for similar short time- 
scale flaring activity during this outburst. We examined the evo- 
lution of the SPI and ISGRI individual detector count rate on 
time scales of a few seconds (and down to 0.1 s with ISGRI) 
around the peak of the outburst but found no evidence for the 
presence of such short flares in our data. PICsIT spectral timing 
data were also checked for the presence of very short flares. With 
the 16 ms time resolution nothing was found in 8 energy bands 
between 208 keV and 2.6 MeV. 

In the context of previous observations, the 2006 September 
event is therefore not extreme in terms of luminosity. The bolo- 
metric flux for our e-folded power law model using the SPI 
data is 8.5 x 10~ 8 and 1.2 x 10~ 7 erg cirT 2 s -1 for the aver- 
age of rev. 482 and its peak, respectively. This corresponds to 
(2.8-3.9)%(10M o /M)L E . These Eddington ratios are still be- 
low those seen in the soft state of Cyg X-l, ^ 5%(10M Q /M)L E 
(Z02), unlike the case of the shorter flares, which L substantially 
exceeded this value. 

However, the INTEGRAL flare is exceptional in terms of du- 
ration. The source stayed in this bright state for several days, and 
the 10-hour ouburst is among the longest recorded. This suggests 
a relative stability of this bright state. The fact that this flare was 
below the flux level typical of the soft state could be related to its 
long duration. Were it above the soft-state level, it would, most 
likely, have triggered a state transtion to the soft state during the 
2-day time scale, perhaps like the one described in Malzac et al. 
(2006). 

We note that another flare, indicated in Fig. [T]by the blue 
arrows, was recently reported by INTEGRAL (Neronov et al. 
2008). It was, however, substantially weaker (it did not even cor- 
respond to a local flux maximum of the ASM/BAT light curves), 
and with a shorter time scale than the 2006 event. On the other 
hand, an event with the amplitude approximately equal to the 
2006 one, but still substantially shorter (of a ~ 1-d duration), took 
place on MJD 54690 (2008 August 12), as shown by the green 
arrows on Fig. Q] It appears to correspond to the maximum of 
the superorbital cycle of Cyg X-l, currently with the period of 
~300 d (Rico 2008). It occured, however, 35 d later than the date 
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E [keV] 

Fig. 8. Comparison of the INTEGRAL spectrum corresponding 
to the peak of the flare (the SPI data are shown in blue, ISGRI in 
cyan and PICsIT in red crosses; the data have been rebined for 
clarity), the average Cyg X-l hard-state spectrum (McConnell et 
al. 2002) from the CGflO/OSSE/BATSE/COMPTEL including a 
BeppoSAX spectrum at low energies (the green histogram fitted 
by the solid curve), and the brightest previously recorded ks- 
time scale flare (the dashed 15-300 keV curve). For that last 
spectrum, we show the fit to the average over the period C of the 
flare 990421A (table 2 in Golenetskii et al. 2003). Note that the 
peak flux of that flare was still 1.8 times higher than that of the 
average spectrum shown here. For clarity the INTEGRAL data 
have been rebined and normalised to the SPI spectrum according 
to the best fit e-folded power-law model. The ISGRI and PICsIT 
fluxes were divided by 0.84 and 0.85 respectively (the PICsIT 
normalisation factor is poorly constrained). 



predicted by Rico (2008), indicating a quasiperiodic character of 
the superorbital modulation. 

Several other sources like XTE J 1739-302, Vela X-l and 
IE 1145-6141, which also exihibit relatively short time-scale 
(hours) flares, share with Cyg X-l the fact of being bright X- 
ray binairies with a giant companion where accretion proceeds 
trough the wind of the companion star (Negueruela et al. 2006; 
Smith et al. 2006). At least some of these bright events are likely 
to be caused by an increased mass accretion rate triggered by the 
variability of the wind. 

The fact that the spectrum observed in this luminous state is 
very similar to the hard state spectrum suggests that the emission 
mechanism is similar. In the hard state, those spectra are usually 
believed to form close to the black hole in some kind of hot 
geometrically thick, optically thin accretion flow. Although the 
simplest advection dominated accretion flow solutions (Narayan 
& Yi 1995) may not hold at such high luminosities, other forms 



8 



J. Malzac et al.: An intense state of hard X-ray emission of Cyg X-l 



of luminous hot accretion flow solutions are known to exist and 
to be sufficiently stable (Yuan 2001, 2004; Yuan et al. 2007). 

Alternatively, the emission could be produced in an outflow- 
ing, optically-thin, patchy corona atop a standard geometrically- 
thin, optically-thick accretion disc (Beloborodov 1999; Malzac 
et al. 2001). The possibility that the hard X-ray emission during 
the flare arises in a jet due to a microblazar activity as suggested 
by Romero et al. (2002) is ruled out as this model predicts a 
softer spectrum during bright flares. Also, the flare periodicity 
predicted by that model has been shown to be absent in the data 
(Lachowicz et al. 2006). 

The INTEGRAL spectrum at the time of the MAGIC detec- 
tion was a typical hard state spectrum, and similar to the peak 
spectrum (during the non-detection). In fact, we do not find any 
unusual features in the X-ray light curves and spectrum around 
the time of the MAGIC detection, except for the high X-ray lu- 
minosity that lasted for a few days. 

As discussed in Section [3] we do not find conclusive evi- 
dence for a non-thermal high energy component in the present 
data. However, as already mentioned, deeper exposure observa- 
tions performed when the source was in the hard state have de- 
tected such a hard tail. It is usually attributed to the presence of a 
small fraction of non-thermal electrons in the hot comptonising 
medium. In the case of Cyg X- 1 , this non-thermal component ex- 
tends at least up to a few MeV (McConnell et al. 2002, see Fig. 
|8}. The TeV emission reported by A07 could in principle be re- 
lated to such non-thermal coronal electrons. If we extrapolate the 
reported TeV emission down to the MeV range assuming T = 2.5 
(a characteristic non-thermal index in Cyg X-l, e.g., Gierlinski 
et al. 1999), we match the SPI residuals at -500 keV (Figs.|6H7]>- 
This spectral index is very close to the range of F = 3.2 + 0.6 
measured by the MAGIC detector (A07). However, it is highly 
unlikely that the spectrum from < 1 MeV up to 1 TeV is a sin- 
gle power law. Even if the intrinsic spectrum is a power law, 
reprocessing by e* pair production events will change its shape. 
The issue of pair absorption and pair reprocessing is, however, 
beyond the scope of this paper and calculations of these effects 
will be presented elsewhere. 

Alternatively, the TeV emission could arise in shocks located 
in the region where the outflow originating close to the black 
hole interacts with the wind of the star. Such an outflow is in- 
deed observed in Cyg X-l in the form of the compact radio jet. 
Other forms of outflows such as an accretion disc wind could 
also be responsible for this emission. This outflow colliding with 
the star wind may power strong shocks accelerating the particles 
responsible for the TeV emission. In this scenario an increase of 
the wind strength would increase the mass accretion rate onto 
the black hole, leading to the enhanced hard X-ray emission. 
Simultaneously the stronger wind interacts more efficiently with 
the (possibly stronger) outflow, making the TeV emission de- 
tectable. Still, pair absorption in the photon field of the compan- 
ion (e.g., Dubus 2006) remains to be a major issue. 
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